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Abstract 

In this contribution, a design of a synthetic cahbration genetic circuit to char- 
acterize the relative strength of different sensing promoters is proposed and 
its specifications and performance are analyzed via an effective mathematical 
model. Our calibrator device possesses certain novel and useful features like 
modularity (and thus the possibility of being used in many different biolog- 
ical contexts), simplicity, being based on a single cell, high sensitivity and 
fast response. To uncover the critical model parameters and the correspond- 
ing parameter domain at which the calibrator performance will be optimal, 
a sensitivity analysis of the model parameters was carried out over a given 
range of sensing protein concentrations (acting as input). Our analysis sug- 
gests that the half saturation constants for repression, sensing and difference 
in binding cooperativity (Hill coefficients) for repression are the key to the 
performance of the proposed device. They furthermore are determinant for 
the sensing speed of the device, showing that it is possible to produce de- 
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tectable differences in the repression protein concentrations and in turn in 
the corresponding ffuorescence in less than two hours. This analysis paves 
the way for the design, experimental construction and validation of a new 
family of functional genetic circuits for the purpose of calibrating promoters. 

Keywords: synthetic genetic circuits, synthetic biology, calibration, gene 
promoter, effective modeling of gene circuits, parameter analysis 



1. Introduction 

One of the fundamental principles of synthetic biology is the construction 
of biological standardized parts and devices which are interchangeables. A 
proper characterization of these parts and devices appears as a key issue in 
order to make them reusable in a predictive way. In the recent past scien- 
tists have witnessed several initiatives towards the design and fabrication of 
synthetic biological components and systems as a promising way to explore, 
understand and obtain beneficial applications from nature. For instance, in 
the post genomic era one of the most fascinating challenges scientists are 
facing is to understand how the phenotypic behaviour of living cells arise out 
of the properties of their complex network of signalling proteins. While the 
interacting biomolecules perform many essential functions in these systems, 
the underlying design principles behind the functioning of such intracellu- 
lar networks still remain poorly understood js], |l3[ . Several initiatives have 
been reported in this line of thought to uncover some key working principles 
of such genetic regulatory networks via quantitative analysis of some rel- 
atively simple, experimentally well characterized, artificial genetic circuits. 
It has been shown that custom made gene-regulatory circuits with any de- 
sired property can be constructed from simple regulatory elements 0] . These 
properties include bistability, multistability or oscillatiory behaviour of ge- 
netic circuits in various microorganisms such as bacteriophage switch jsf or 
the cyanobacterium circadian oscillator [g!] . As one example, the genetic tog- 
gle switch, a synthetic, bi-stable gene-regulatory network in Escherichia coli, 
was shown to provide a simple theory that uncovers the conditions necessary 



for bi-stability |12| . Further, artificial positive feedback loops (PFLs) 
have been used as genetic amplifiers in order to enhance the responses of 
weak promoters and in the creation of eukaryotic gene switches [Mj- Sayut 
et al. demonstrated the construction and directed evolution of two PFLs 
based on the LuxR transcriptional activator and its cognate promoter, Pluxl 
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These circuits may have apphcation in metabohc engineering or gene 
therapy that requires inducible gene expressions j^, . 

The desired performance of these synthetic networks and in turn the 
resultant phenotype is strongly dependent on the expression level of the cor- 
responding genes, which is further controlled by several factors such as pro- 
moter strength, cis- and trans-acting factors, cell growth stage, the expression 
level of various RNA polymerase-associated factors and other gene-level reg- 



ulation characteristics [ll|, |13|. Thus, one important ingredient to elucidate 
gene function and genetic control on phenotype would be to have access to 
well-characterized promoter libraries. These promoter libraries would be in 
turn useful for the design and construction of novel biological systems. There 
have been several initiatives to control gene expression through the creation 
of promoter libraries [21,0] ■ Alper et al., [l[ have reported a methodology to 
develop a completely characterized, homogeneous, broad-range, functional 
promoter library with the demonstration of its applicability to analysis of 
genetic control. 



Since Miller published [16| a proposal for a measurement standard for 



/3-galactosidase assays, yet much work has been done with no conclusive 



standard being established |17l. Il8l. Il9[|. The main goal in calibration is mea- 



suring a query value up to an established standard. A good device should 
be unique, reliable and easy to use; additionally it should circumvent, to 
all possible extent, any noise that could alter the measurement. Recently a 
methodology 20(] has been reported to characterize the activity of promoters 
by using two different cell strains. In the present study we propose the use 
of a synthetic gene regulatory network as a framework to characterize dif- 
ferent promoter specifications by using a single-cell strategy. In this context 
characterization stands for evaluating the parameters of a query promoter as 
compared to a standard promoter acting as a scale. The proposed device, 
the promoter calibrator, works on the principle of comparing a specific input 
signal which will be sensed by promoters of different sensing strengths and, 
as an output, produces fluorescence of specific colours which allows quanti- 
fying the relative strength of the promoters. Analyses were carried out in 
order to find out relevant model parameters and the corresponding range of 
model parameter values which are compatible with the performance of this 
calibrating biological design over a spectrum of given input . 

This contribution is organized as follows: in the first part, "Design", the 
structure and working principle are explained and the mathematical model 
resulting from the construction is established. In section [31 "Numerical Anal- 
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Figure 1; Design of the proposed promoter calibrator. It is composed of two promoters 
(with two parts each: a sensing and a repressed domain) one of the sensing promoters 
is the device promoter and the other is the query promoter. The repressed domains are 
controled by the two repressors proteins (x and y). Each promoter is inhibited by the 
repressor which is transcribed from the opposing promoter. Fluorescence proteins levels 
will be proportional to repressor protein levels, which, in turn, will be promoted by the 
sensing promoters. 



ysis of the System" , we analyze the dynamics of the model equations in regard 
to its stability, functional parameter regions and sensitivity or robustness vs. 
the change in certain key parameter values. In the following section, a proof 
of concept design is proposed in order to choose the right parameters to actu- 
ally perform the experimental validation of our concepts and have a system 
that gives a clear and stable signal that can be interpreted. Finally, the 
conclusions resulting from our paper are exposed. 

2. Design 

2.1. Biological principles 

Our promoter calibrator is composed of two promoters (each with two 
parts: a sensing and a repressed domain), two repressors proteins and two 
fluorescent protein outputs (see Fig. [1]). Each promoter is inhibited by 
the repressor, transcription of which is promoted by the opposing promoter. 
Fluorescence protein levels will be directly related to repressor protein lev- 
els, activated in turn by their sensing promoters. Hence, different sensing 
strengths will cause a difference in the expression of the fluorescence pro- 
teins, detectable by means of single cell fluorescence as changes in the color 
patterns of the individual cell or cell sample. 

In our scheme, one of the sensing promoters acts as the device promoter 
to which the strength of a given query promoter is quantitatively compared. 
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The main use of this device is to characterize different promoter specifications 
(sensing affinities and cooperativities) compared to some standard. One of 
the main usefulness of this design hes in the potential modularity of the sys- 
tem: by changing the sensing part of the promoters, other sensing promoters 
could be calibrated; this change can be carried out by a simple, straight- 
forward cloning step. Modularity also boasts the potential of this device as 
it can be implemented in a potentially unlimited set of systems. 

2.2. Mathematical model 

The behaviour of the proposed promoter calibrator can be understood 
via an effective mathematical model. The model is considered to be effective 
as transcription and translation have been modeled as a lumped reaction. 
The separation of transcription and translation otherwise involves a response 
delay. We seek to classify dynamic behaviors depending upon the change in 
model parameters and determine which experimental parameters should be 
fine-tuned in order to obtain a satisfactory performance of our device. 

The time dependent changes in repressor and sensing protein (input) con- 
centrations is shown in equations ([T]|3]). Subsequent to the biological design, 
reporter protein concentrations are directly related to repressor protein con- 
centrations. 



dx {tj 1 

«1 / xm / - PxX + -fx, (1) 



dt 



1 + + 



n2 

Ps 



dy Ifca, ^ ^ . . 

= "2 ^^n^-PyV + ly, (2) 



l+(f 



dps 
dt 



-PpsPs- (3) 



The device and query promoters activate the production of repressor pro- 
tein X and y, respectively, and their concentration is related directly to the 
concentration of fluorescence proteins. Thus these variables will be treated as 
equivalent from the modelling point of view. Parameters ai and 02 represent 
the effective rate of synthesis of repressor proteins x and y, respectively; a is 
a lumped parameter that takes into account the net effect of various activ- 
ities such as RNA polymerase binding, RNA elongation and termination of 
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transcript, ribosome binding and polypeptide elongation and will be modi- 
fied by repression and sensing effects. The (3^, (3y and (3p^ arc the degradation 
constants of repressor protein x, repressor protein y and sensing protein ps, 
respectively. The sensing protein concentration ps will depend on the sensed 
input, will be easy to change in a given experiment and is used as the main 
input variable in our calibrator experiments. It is important to note that a 
slow rate of degradation is assumed for the sensing protein, implying a nearly 
constant level over a reasonable experimental time interval. Basal level rates 
of synthesis of proteins x and y are denoted by and jy, respectively. 

Repressor and sensing responses are assumed to follow Hill equation dy- 
namics: promoter-binding monomers form multimers by positive allosterism 
and attach to its cognate promoter with saturating behaviour. Binding coop- 
erativities are described by Hill coefficients fix and Uy for repressor domains 
corresponding to x and y respectively, and rii and n2 for sensing domains 
corresponding to device and query promoter respectively. The extent of the 
saturation rate is described by half saturation constants or Michaehs con- 
stants, denoted by parameter and ky for repressor domains corresponding 
to X and y respectively and ki and k2 for sensing domains corresponding to 
device and query promoter respectively. The total number of promoter sites 
is assumed to be conserved and the total concentration of both promoters is 
chosen to be identical. 

In our construction, the crossrepressing part will be kept unchanged while 
different sensing domains may be attached to it. The aim is to establish 
a protocol to accurately quantify differences between the sensing promoter 
parameters (q;i,2, ^1,2)- Crossrepression parameters {kx,y, Px,y and nx,y) are 
structural parameters that must be chosen in such a way that the fluorescence 
response of the system gives us stable, sensitive and robust indication about 
the quantitative relations between the sensing promoter parameters. The 
dynamic analysis of the system will help us to take the right decisions on 
which are the most appropriate values for these structural parameters. The 
next sections are devoted to the dynamical analysis in order to determine the 
sensitivity and robustness of the system for different ranges of the structural 
parameters. 

The commercial software package Mathematica (Wolfram), was used for 
model development and simulation. In the numerical calculations we have 
used the following dimensionless variables: 
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therefore, the units in the plots of the figures in this work are given in units 
of kx or ky for the x and y repressor proteins concentrations and time in units 
of For the adimensional variables, Equations (Il]l2]) take the form: 



\ "2 



where R is the ratio 4^. 

Px 



3. Numerical analysis of the system 

The simplifying assumption of considering sensing proteins for which the 
degradation constant f3p^ is much smaller than the rest (/3p^ <C (3x, (3y) was 
made in order to classify the possible dynamic scenarios of our model. Given 
this assumption, in a first order of approximation we have. 



^ = -^.Ps^O. (11) 

In such approach, the concentration of sensing protein ps is constant 
during the evolution time of the rest of the internal variables of the system. 
This assumption leads to a system of two autonomous coupled non-linear 



7 



ordinary differential equations dependent on tlie variables x and eqs. (P- 
[TU]) . in which Ps is fixed although it can be easily changed within a given 
experiment. This is not true for the rest of parameters which are more 
difficult to modify in a given experiment. This approximation transforms 
the system into: 



dY 1 

— = a'^- —-RY + %. (13) 

where the new parameters a[ (effective transcription factors) are given by 
the following expression: 

«' = ^Yvkf' ^^^^ 

In the limit in which the constants k^^y, lix,yi lx,y are equal, this equations 
describe the biological equivalent of an electronic comparator^ that is, a device 
which compares two voltages or currents and switches its output to the larger 
signal. In the biological equivalent, our comparator would select for the larger 
of the two a's, as exemplified in Fig. [21 which represent the evolution of the 
system for the cases in which the query promoter has a higher and lower 
effective transcription factor compared to the device promoter, respectively. 

In any case, our aim is to construct a device, termed a calibrator, which 
not only selects the stronger affinity but also allows quantifying the rela- 
tive strength of both promoters. Although the comparator is a fundamental 
part of this device, a deeper understanding of the dynamics of the system is 
required for its application as a calibrator device in real biological environ- 
ments. 

3.1. Dynamic analysis of the calibrator 

The dynamical analysis of the system given by Eqs. f lT2lfT3|) requires the 
determination of its steady state solutions and their linear stability. The 
steady states {xss.yss) are given by the intersection of the null clines: 



Fi{X,Y) = a;^^^-X + 7. = 0, (15) 
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Figure 2: Typical response of the proposed promoter calibrator. In the upper figure the 
concentration of the x protein (solid line) in the steady state is higher while in the figure 
below the concentration of the y protein (dashed line) is higher. 



and 



F^{X,Y) = 4-_L__i?y + ^^ = o. (16) 

The analytical solution of Eqs. ( |T5lfT6l) cannot be obtained, hence nu- 
merical methods must be used. The linear stability of the steady states is 
determined by the sign of the eigenvalues of the Jacobian matrix, 

M = ( m ii] (17) 

\ ax dY / X=Xss,Y=Yss 



which are given by 



1 + R . 1 



A± = -^±-v/(i?-l)^ + 4A, (18) 



From the analysis of the previous equations fll5H16p . we deduce that, for 
the positive steady state solutions {Xgs > and Yss > 0), the following 
mathematical constraints hold: a[ > Xgs — 7x > and ^2 > YssR — 7y > 0, 
respectively. Thus, taking into account fll8til9p . we observe that A > and 
A_ is always negative. However, A+ can be either negative, for A > R, 
or positive, for A < R, resulting in either stable nodes (sinks) or unstable 
saddles, respectively. The condition A = i? is satisfied at certain critical 
values of the parameters at which precisely one of the steady state solutions 
of the system changes its stability. 

In order to highlight the specific aspects of the calibrator dynamics, we 
will in the following sections consider a number of special cases. Specifically 
we will examine the (fully) symmetrical calibrator, a[ = 02 = = Uy = 

n, kx = ky = k, f3x = (3y ^ R = I and = % = 7, and the partially 
symmetrical calibrator, with the same specifications except that a'l and a'2 
may differ. At the end of the section some general considerations about 
dynamics of the system in the most general case will made. 

3.2. The fully symmetrical calibrator (a[ = cig = «V 

From the analysis of Eqs. ( lT5]fT6l) it is shown that there is always a fixed 
point with Yss = Xss and that there exists a minimum value of Xm such that 
for parameters resulting in Xss > X^, three steady states exist, otherwise 
only one. 

Using a' as free parameter and taking fixed values for the rest, i.e., n, R 
and 7, the condition A = i? = 1, together with Eq. (fT5|) . allows to obtain the 
critical values and Xm that characterize the appearance of the bifurcation, 
namely: 



. _ - XJ\^ - X^ + a'J^ 

whose values can be obtained by numerical methods. For example, for n = 2, 
k = 80, (3 = 0.069 and 7 = 0.1, yields a'^ = 11.24 and = 81.46 or, in 
the dimensionless variables: X = 1.018 and a' = 2.036. Figure 121 shows the 
bifurcation diagram for Xss as function of a showing that for a > am there 
are three steady states. 

This analysis shows that the (fully) symmetrical calibrator possesses three 
fixed points for d[ > a'^: a saddle [xli) with Xgs = Ygs, and two sinks, one 
with Xss > Yss and another one with Xgs < Yss, referred to as x/j and x^. 
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respectively. This behaviour is typical of the occurrence of a (supercritical) 
pitchfork bifurcation and bistable behaviour. 

Regarding the possible trajectories of the dynamic variables, Figure H] 
illustrates the phase plane of Eqs. ( !T2lfT3|) . where the steady states are located 
at the intersection of the null clines eqs. fll5til6p represented by dashed lines. 
The solid lines are the stable (W^) and unstable (W^) manifolds of the 
saddle fixed point xli. The stable manifold divides the phase plane 
in two regions, the first and second octants corresponding to the attraction 
basins of the sinks xr and xl, respectively. Different possible trajectories in 
the phase plane are depicted for a given number of initial conditions, where 
the arrows indicate the flow direction. 

In a calibrator experiment the initial value of the repressor protein con- 
centrations X and y would be zero and hence the phase plane trajectories 
would depart from the origin in Figure 4. For values of a' larger than a'^, 
the system becomes unpredictable, as small perturbations in the trajectories 
would potentially push the system into any of the attraction basins of the 
sinks xTr and xi. 

3.3. The partially symmetrical calibrator 

We consider now the more general scenario in which a[ and 02 may differ 
being the rest of variables equal (n^ = Uy = n, = ky = k, (3^ = f3y = (3 
and 7x = % = 7)- The condition A = R which characterizes the occurrence 
of the pitchfork bifurcations now reads: 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 4: Phase plane, showing the unstable equilibrium point (the point where the two 
dashed lines touch in the center) and the two steady state solutions (points where the 
dashed lines touch close to each axis). The arrows show the path the system would do 
starting from any point in the phase space. 



X Y a'a' ^ ^ 

^^ss^ ss^x y 

that shall be solved together with Eqs. fll5til6p for the fixed points of the 
system. 

Fig. O shows the result of the numerical simulation of the resulting system 
of equations (with initial conditions X = F = 0) by slightly changing the 
value of with respect to The figure shows the results of different 
simulations for a'l =3.0, n = 3 and = ea[ with e =0.5, 0.6, 0.7, 1.0, 
1.5. The results for e < 1 are the points in the right down corner of the plot. 
One can see that these points positions are very insensitive to the value of 
There is only one point in the center of the plot, which corresponds to 
a[ = 0(2, it is the unstable saddle, and small perturbations in the system will 
drive the system away from this solution to either of the other two steady 
state solutions. Once e > 1, the system goes to the solutions where Ygs > Xss 
which are represented by the points in the upper left corner. For these points 
the maximum value of a' is growing and one can observe that the solution 
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Figure 5: Results for the simulation of the partially symmetric calibrator. On the right, 
close to the x-axis {a'2 < a'l for these points), there are many points at the same posi- 
tion, showing that, for the parameter region where a bifurcation happens, the solution is 
insensitive to the value of the weakest between the two a's. 



is sensitive to this value. So the steady state solution into which the system 
falls is only sensitive to the bigger value between a[ and ai^ and changes in 
the smaller among these two parameters has no sensible effect in the final 
solution. 

For the case in which the calibrator falls within the region of bistability, 
if a'2 < a'l the orbits departing from the origin of Fig. |5] would fall within the 
attraction basin of solution xr. It is nevertheless observed that x"r is quite 
insensitive to the actual (y'2/a'i ratio. In consequence, the system would show 
a stable but rather insensitive response to different query promoters. On the 
other hand, if a[ < a'2-, the orbits departing from the origin would fall within 
the attraction basin of solution xl, which changes appreciably as a function 
of the a'2/ a'l ratio. Thus the system would not only be stable, but also rather 
sensitive to changes in the effective query promoter affinity. It should be kept 
in mind that the sensing protein concentration, p^? can be used to modify 
a'l, a'21 which changes from unity to 0.1^2 as Ps changes from zero to infinity 
and therefore the ratio a'2/a'i changes with p^. 

We can also define the fluorescence ratio as the ratio of X/Y if X < y 
and Y/ X iiY > X. This will be the intensity ratio of the two fluorescences 
once the system reaches stability. In Fig. |6] we show a plot of this ratio for 



13 



different values of a'2lo'x- This ratio grows until it reaches its maximum when 
and then it decreases. Another observation about this parameter 
is that the bigger a'^ is, the less sensible to the ratio a.^ja'^ the fluorescence 
ratio will be. 

3.J^. The calibrator dynamics in the general case 



The theorem of Andronov and Pontryagin [21| states that Eqs. 
in the symmetrical case are structurally stable, since every fixed point is hy- 
perbolic (its eigenvalues have a non-null real part) and there are no orbits 
connecting two saddles (since there is only one). Structural stability im- 
plies that the phase plane topology is preserved under small perturbations 
of the parameters. Hence, the phase plane of Eqs. fll2m3p in the case that 
a'^ ~ a'y, Tlx ~ Uy, kx ~ ky, (3x ~ f3y and •jx ~ is topologically equiv- 
alent to that shown in Fig. 4, meaning that there is a continuous function 
(homeomorphism) between both phase planes. 

Changing the ratio of other structural parameters of the calibrator has 
similar results as in the partially symmetrical case. For a given range close 
to the value 1 for the ratio of each parameter ratio {ux/y, Px/y, •••) the 
bifurcation appears while far from the value 1 the bifurcation cannot be 
seen. The range is usually bigger, the bigger the values for a[ 2 ^^e. In Fig. [7] 
we show, as an example, the range where the bifurcation appears for different 
values of (3x/ (3y. 

If i? < 1 the orbits departing from the origin (X = Y = 0) would fall 
within the attraction basin of solution xl, on the other hand if i? > 1 the 
orbits departing from the origin would fall within the attraction basin of 
solution xr. 

3.5. Calibrator performance analysis: robustness and response time 

In order to use this system to measure the relative strength between two 
promoters, one should keep in mind two factors. The first important factor 
is the right choice for the parameters of the repressor proteins and device 
promoter in order to have a robust system, that gives a stable response that 
can be easily interpreted. Second, is the time response of the device, that 
means, how long does the system needs to reach its steady state solution. 

When the equations are written in the dimensionless form, the parameters 
kx and ky do not appear explicitly, see eqs. fl9l lT0|) . These parameters appear 
implicit in the definition of the variables X and Y and in the 7 parameters 
(which have small influence in the dynamics of the system). By choosing 
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Figure 6: Upper plot: Fluorescence ratio for different values of oi'y/a.'^ (a^ =2.5). The blue 
points are solutions where X >Y and in the red points Y > X. Lower plot: Fluorescence 
ratio for different values of a'y/a'^ and for different values of a'^ (Solid line:Q:^ =2, dashed 
line:a^ =3, dotted line:a^ =4). 
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Figure 7: (Color online) Position of Xss for different values of /3y//3x- In the black curve 
a'^ — a' —2, in the blue curve a'^ — a' —3 and in the red one a'^ = a' =4. 



kx = ky the results will be easier to interpret since the fluorescence is directly 
related to the concentrations of the proteins x and y and, by setting k^ = 
ky, the fluorescence intensity ratio [X/Y and Y/X) and the fluorescence 
intensity difference — F |) will be directly proportional to these parameter 
calculated with the real protein concentrations. 

An experiment made with the calibrator would consist of cloning a plas- 
mid with the calibrator genetic circuit assembled with the device promoter 
(whose parameters one have to choose) among known ones and with a query 
promoter whose parameters are unknown. The plasmid should be inserted 
in cells in solutions of the signaling protein at different concentrations ps- 
Each promoter is modeled through two parameters, 0:1/2 and ki/2, 1/2 stand 
for device/query promoter. While at low ps concentrations both promoters 
are weak and give a weak fluorescence response, at high ps concentrations, 
both promoters are saturated and their strength is maximal. From the flu- 
orescence intensities at these high concentrations of the signaling protein it 
is possible to establish the relative strength of the two promoters a2/ai- In 
figures [8] and [9] we show plots of the fluorescence difference defined as |X — F| 
and the fluorescence ratio X/Y for three different values of di and varying 
a2 at high signaling protein concentrations (the effective strength of both 
promoters is maximum). 
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Figure 8: (Color online) The fluorescence difference for different values of ai as a function 
of a2- Note that for values of a2 sufficiently higher than ai the fluorescence difference 
increases linearly with the value of 6:2 ■ 



The first thing to note from figures [8] and |9] is that, if the query pro- 
moter is stronger than the device one, the device fluorescence (X) will be 
strongly suppressed, and the fluorescence intensity coming from the query 
promoter is proportional to its strength (the response of the system is lin- 
ear). That means, choosing a weak device promoter, one can establish the 
relative strength of other promoters by a simple proportionality law given by 
the linear response plotted in figure El 

At each different Ps concentration, the effective strength of the device and 
query promoters is different, see eq. ( I14p . The parameter that distinguishes 
two promoters, with respect to the ps concentration, is their Michaelis con- 
stants, ki^2- The parameters ki^2 mark the rhythm at which the effective 
strength of each promoter grows. If a promoter has a small value of k, at low 
Ps concentrations of the signaling protein, the promoter is already acting at 
full strength, while for high values of k the promoter saturates only at high 
values of ps- We have already established to choose a small value for the 
device promoter ai, so we expect the query promoters to have 0:2 > ai- If 
/c2 < ki, the effective strength of the query promoter is always bigger than 
the relative strength of the device one, and in the experiment one observes 
that the luminosity associated with the query promoter is stronger for any 
value of the signaling protein concentration p^. On the other hand, if one 
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Figure 9: (Color online) Upper plot: the fluorescence ratio X/Y as a function of a2/ai 
for different values of ai. One can clearly see that for similar values of ai and a2, when 
the bifurcation occurs, the system goes to a state where the repressor protein of the 
stronger promoter completely dominates the system. Lower plot: Detail of the region 
where 0:2 > ai. 
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chooses a small value for ki , already at low Ps concentrations the strength of 
the device promoter saturates, and if ki is small enough it saturates before 
the effective strength of the query promoter reaches a value bigger than ai. 
In this situation one would observe at low concentrations of ps the lumi- 
nosity of the device promoter stronger than the one coming from the query 
promoter. Then, at some critical value of Ps = Psc both strength are equal 
and for p^ > Ps^ the stronger fluorescence is the one from the query promoter. 
For ni = n2, the value of ^2 given in units of ki as a function of psc (also in 
units of ki) is given by: 



In figure [TO] we show a few examples of results one might expect for 
different values of /c2- 

So, the construction of the calibrator device, as we present it, would be 
the following: first one chooses a very weak promoter which has a small 
Michaelis constant to act as the device promoter in the calibrator. Second 
step is to define a standard, to choose a known promoter, clone the calibrator 
device with it as query promoter and perform a measurement of the fluores- 
cence intensity of this standard promoter at high ps concentrations. This 
fluorescence intensity is the standard one, to which we can compare other 
promoters. Now performing the experiment with another promoter acting 
as query promoter one obtains another value for the luminosity that we can 
compare with the standard one. The higher or lower this luminosity is with 
respect to the standard, the stronger or weaker the promoter is compared 
with the standard, so one can establish the value of 02 • Knowing 02 one can 
perform the same measurement for different Ps concentrations in order to 
establish the critical value of ps where the query fluorescence becomes higher 
than the device one. Knowing the value of psc it is possible to establish the 
value of k2 by means of eq. [22] (assuming both promoters have the same n). 

Now that we have established the ideal parameters for the device promoter 
(weak strength and small Michaelis constant) and set kx = ky and (3^ = Py 
the last important factor is the time response of the system. 



Psc 




(22) 



(23) 
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In figure [TT] we show plots for tlie tf, the time tlie systems needs to 
reacli its steady stat^ for different values of ai. One observes tliat tlie time 
response of the system has a peak with the maximum around 30/3"^ when 
the effective strength of both promoters is equal and then it goes to a rather 
stable value close to 7/3~^. For a realistic value of (3x like 0.069 min~^ the peak 
value for t/ is 7 hours, while for most of the measurements (the calibrator at 
different Ps concentrations) this time should be around two hours. 

4. Conclusions 

In the present study we have proposed a biological device that works as a 
promoter calibrator in which the strength of a collection of query promoters 
can be measured against the strength of a device promoter. Some of the key 
features of the proposed design are its single cell character, high modular- 
ity and handy construction: a unique molecular cloning permits the change 
of the promoter ready to be calibrated. The designed performance of the 
proposed biological device has been demonstrated by means of an effective 
mathematical model. The sensitivity analysis of the model shows that there 
is a sensible relation between the relative promoter strengths and the final 
steady fluorescences measured by the system. 

Furthermore, a response time analysis shows that the device can produce 
a large difference in the repression protein concentrations and in turn in the 
corresponding fluorescence in approximately two hours. 

Finally our promoter calibrator principle may lead to an improvement in 
the modeling and characterizations of systems in Synthetic Biology, which 
frequently rely on arbitrarily characterized, or even non-characterized, pro- 
moters. 
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Figure 10: In all plots ai- 
and in the bottom plot k2 = 



-ki—1 and a2—2. In the upper plot k2—2, in the center fc2=3 
=4. The values for psc are respectively: V2, V7 and -^14. 
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Figure 11: Response time of the system for different values of a\ and a.2- For values of 
the a's that the system presents a bifurcation, the response time can be large because the 
system spend time in its non-equilibrium solution. 
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